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Time Reversal Communication in Multi-Path Fading Channels with Pinholes 
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The paper presents an analysis of the time reversal in multi-path Rayleigh-fading channels with 
N inputs (transmitters) and M outputs (receivers) . The main issues addressed are the condition of 
statistical stability, the rate of information transfer and the effect of pinholes. The stability condition 
is proved to be MC <C N e gB for broadband channels and M <C N e g for narrowband channels where 
C is the symbol rate, B is the bandwidth and N e s is the effective number of transmitters. It is 
shown that when the number of layers, n — 1, is relatively low compared to the logarithm of numbers 
of pinholes N e g is given by n~ times the harmonic mean of the number of transmitters and the 
numbers of pinholes at all layers. On the other hand, when the number of layers is relatively large the 
effective number of pinholes diminishes exponentially. The energy efficiency is shown to be optimal 
when the power supply is set to the noise level times BN e h and that the maximal information rate 
is roughly BN c g when the stability condition is violated. 

PACS numbers: 46.65. +g, 43.20.+g, 42.68.Ay 



I. INTRODUCTION 

Time reversal (TR) of waves has received great at- 
tention in recent years and been extensively studied 
for electromagnet ic Bl , [2fj| . |33 | as well as acoustic 
propagation (see |14| and the references therein). A 
striking effect of time reversal in randomly inhomoge- 
neous media is the superresolution of refocal signals 
1 ^3 which implies low probability of intercept and 
holds high potential in techn olog ical applications such 
as communications 0, [lU, |10| - [l7| - 

An issue prior to superresolution, however, is the 
condition of stability, namely how many antennas and 
how much bandwidth one needs to achieve statistical 
stability in TR so that the received signals are nearly 
deterministic, independent of the channel statistics? 
In this note I answer this question for multi-path 
Rayleigh fading channels, with multiple inputs and 
multiple outputs (MIMO), commonly used in wireless 
communication literature, see, e.g. 23] and I analyze 
the effect of multi-layer pinholes in Section Ivl 

In the MIMO-TR communication scheme||, 0, 
the M well-separated receivers first send a pilot sig- 
nal to the TV-element TRA which then uses the time- 
reversed version of the received signals to modulate 
the data symbols and retransmit them back to the re- 
ceivers. One of the main results obtained here is that 
the time reversal process is statistically stable when 

MC <C N c gB, for broadband channels (1) 
M <C N c g, for narrowband channels (2) 

where C{< 2B) is the symbol rate, B is the bandwidth 
and N e s is the effective number of transmitters. In the 
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presence of (n — l)-layer pinholes, I show that the ef- 
fective number of transmitters is asymptotically equal 
to n times the harmonic mean of the number of 
transmitters and the numbers of pinholes of all layers 
when these numbers are greater than 2™ . That is, the 
multi-layer pinholes have a screening effect on time- 
reversal transmission, reducing the effective number 
of the time-reversal elements by at least a factor of 
n. When the last condition is violated, the situation 
is even worse: N c g diminishes exponentially fast as 
the number of layers increases, resulting in a tough 
environment to perform time reversal of good quality. 



The LHS of is the number of degrees of freedom 
per unit time in the constellation of intended data- 
streams while the RHS of Q is roughly the number 
of degrees of freedom per unit time in the channel 
state information (CSI) received by TRA from the 
pilot signals. The latter has to be larger than the 
former in order to reverse the random scrambling by 
the channel and achieve deterministic outputs. The 
stability condition N 1 for narrow-band channels 
or B ^> p c (the coherence bandwidth) for broadband 
channels, when M is small and the pinholes are ab- 
sent, have been previously discussed in 0, [H, 



In Section I take into account the effect of noise 
and analyze the information rate in TR communica- 
tion in noisy channels. I demonstrate a tradeoff be- 
tween stability and spectral efficiency: the maximal 
information rate R ~ BN c ff is achieved when the in- 
equality in is reversed. Also, for a given level 
of noise the energy efficiency is optimized when the 
power level is set to the noise level times N e gB. 
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Let us note that while all the TRA-elements are co- 
ordinated and synchronized the receivers do not have 
the knowledge about the channel and can not coor- 
dinate in decoding the total signals received. As a 
consequence, the co-channel interference from multi- 
ple users can be a serious problem [2^ . An advantage 
of the time reversal scheme is the possibility to use the 
(statistical) stability property to achieve the following 
asymptotic 



FIG. 1: MIMO Channel 



II. TR-MIMO COMMUNICATION 



N 

' " 9(^)^2 Hj k (uj)H* k (uj)duj ~ Sij6(t - n) 

so that Sj(t) ~ Y^iLi m j{ T i)b{t — T i) an d each receiver 
receives the intended symbols without interference. 



First let us review the MIMO-TR communication 
scheme as described in (l2^ . 

The M receivers located at yj,j = 1,...,M first 
send a pilot signal J e luJt g(ui)dLu5(x.j — y^) to the TV- 
element TRA located at — 1,...,N which then 
use the time-reversed version of the received signal 
J e luJt g(oj)H(yj,Xi;uj)duj to encode a stream of sym- 
bols and retransmit them back to the receivers. Here 
H is the transfer function of the propagation channel 
at the frequency 10 from point y to x and g 2 {uj) is 
the power density at uj. Let H(cj) = [Hij(uj)) be the 
transfer matrix between the transmitters and receivers 
where Hij(uj) = H(x.i, y^; to). The reciprocity implies 
that H(w) is symmetric and the relation H*(lu) = 
H(— uj) where * stands for complex conjugation. Let 
us assume that g is a smooth and rapidly decay- 
ing function with effective support of size B such as 

g 2 {uj) = (2 7 r)^ 1 / 2 exp(- l " 2 '^ 12 ). Naturally the rel- 
ative bandwidth B/ujq is less than unity. We have 
chosen the time unit such that the speed of propaga- 
tion is one. 

Let us assume that the separation L between the 
TRA and the receivers is much larger than the spacing 
within the TRA-elements and the receivers. The the 
signal vector S = (Sj) arriving at the receiver with 
delay L + 1 is then given by Q , 0] 

W M „ 

1=1 i=i J 

N 

x£ff Jfc ( w )ff; fc ( w )dw (3) 
fc=i 

where irij(Ti),l = 1,...,W < 00 is a stream of sym- 
bols intended for the j-th receiver transmitted at 
times T; = It,t > 0. In vector notation, we have 
s = Ez=i /e- ,w(f - Tl) .g(w)HHt(w)m(r i )(L where Ht 
is the conjugate transpose of H and 111(77) = (rrij(Ti)). 



III. STATISTICAL STABILITY 

One of the main goals of the present note is to char- 
acterize the stability regime for the important chan- 
nel model of multi-path Rayleigh fading in which Hij 
are i.i.d. C/V(0, 1), the zero-mean, variance-one circu- 
larly symmetric complex-Gaussian random variables. 
For simplicity, I assume that 1 771^(77)! = /i,Vi,L The 
multi-path Rayleigh fading after proper normalization 
is a simplified model for richly scattering environ- 
ment when the spacing within the transmitters and 
receivers is larger than the coherence length £ c of the 
channel. In general, the coherence length is inversely 
proportional to the angular spread |28j | and sometimes 
can be computed explicitly in terms of physical prop- 
erties of the channel 0] • For diffuse waves the coher- 
ence leng th is known to be on the scale of wavelength 

Let us calculate the mean and the variance of the 
signals with respect to the ensemble of the channel. 
We use E to denote the channel ensemble average. 
Consider the quantity HELm. By the Gaussian rule 
for the calculation of moments we have 

EHtfm = iVm (4) 

M 

E|(H(w)H t (w)m) i | 2 = N 2 \m 3 \ 2 + N^\m t \ 2 {h) 

so that ES = NBmJ2Zi e ~ tuJo{t ~ T ' )e ~ B2{t ~ n)2 ■ Let 
t > (2B)^ 1 so that the summation in ES is B- 
uniformly bounded as W — * 00. 

We measure the statistical stability of the signals by 
the normalized variance of the signals at the receivers 

V,(t„) = |E ||^ r Vjirn) = E^fKHES.K)! 2 

Vj, n and say that the signals are stable when 

V i (Tn)-0,Vj ) n. 
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Let P c be the coherence bandwidth of the channel 
such that 

E{(H(a;)Ht(a;)m).(H(a;OHt(a;Om)*} 



N*\ mj \* + NYZi 
0, 



\u-u/\ <p /2, 
\lo - u'\ > f3 c . 



The coherence bandwidth p c is inversely proportional 
to the delay spread and hence the dela y-sp read- 
bandwidth product (DSB) is roughly BP' 1 ]n|, [H , 
|'2.'"{l| . In the diffusion approximation /3 C is given by the 
Thouless frequency DbL~ 2 where Db is the Boltz- 
mann diffusion constant, equal to the energy trans- 
port velocity times the transport mean free path, and 
L the distance of propagation 0] , . 

The broadband, frequency-selective (BBFS) chan- 
nel is naturally defined as having a large DSP, i.e. 
BP' 1 > 1. Since B < uj , w e [uj - B/2,uj a + B/2] 
and — uj are separated by more than p c in a broad- 
band frequency-selective channel. The same holds for 
a narrow-band, frequency-non-selective (NBFN) chan- 
nel defined as B < p c <C ojq. 

Consider the NBFN case first. We have 



M 



Vj (t) « NB^\mi\ 2 \Yl 



e iw T, e -B 2 (t-r,) 2 /2| 2 



;=i 



In view of 10 the stability condition N ^> M for 
NBFN channels then follows easily. The main focus 
of the paper, however, is the BBFS channels for which 
we have instead 



M W 



x — ., v 2 sin %(i — tA 



*-7J 



iui (l' -V)t -B 2 (l-l') 2 r 2 /2 



X X/ 6 ~" V " "'' e 

Several observations are in order. First, due to r > 
(2B)^ 1 the summation over I' in © is convergent as 
W — > oo uniformly in I and B. Second, the summation 
over I is also convergent as W — > oo with the effective 
number of terms ~ BP" 1 . As a result, it suffices to 
consider the case W = 0(BP~ 1 ) or equivalently \t\ — 
i~w\ = 0(P~ 1 ) for which we have the estimate Vj ~ 

NBCj^iLi \ m i\ 2 where C is the number of symbols 
per unit time in each data-stream. It then follows that 



V, 



if and only if NB > MC for BBFS channels. 



The transition to the condition N > M for NBFN 
channels takes place when B ~ C, i.e. r ~ B^ 1 . 

Since iVB is the number of degrees of freedom in 
the channel state information collected at the TRA 
per unit time and MC is the number of degrees of 
freedom in the ensemble of messages per unit time 
the stability condition NB ^> MC can be interpreted 



as saying that in order to recover the deterministic 
messages, independent of the channel ensemble, and 
thus reverse the random scrambling by the channel 
the former must be much larger than the latter. 

A detailed, rigorous analysis of the MIMO-TR 
channel modeled by a stochastic Schrodingcr equa- 
tion, in the paraboli c ap proximation, with a random 
potential is given in |12| . 



IV. RATE OF INFORMATION TRANSFER 

In this section we discuss the information rate 
for a memoryless channel based on the multi- path 
Rayleigh- fading transfer matrix defined above |3l| . 
[l5j. In a memoryless channel an independent real- 
ization of the random transfer matrix is drawn after 
each delay spread. Let us assume as usual that the 
noise at the each receiver is additive-white-Gaussian- 
noise (AWGN) and that the input vector is multivari- 
ate Gaussian and that the channel, the receiver noise 
and the input are independent of one another. 

According to Shannon's theorem the informa- 
tion (in nats) that a symbol can convey on aver- 
age is 2~ 1 ln(l + SINR) where SINR, the signal-to- 
intcrfcrcnce-noise ratio at each receiver, is given by the 
harmonic mean of the SIR, the signal-to-interference 
ratio and SNR, the signal-to-noise ratio, i.e. SINR = 
(SIR -1 + SNR" 1 )" 1 . This formulation assumes that 
the interference is approximately AWGN, like the 
noise at the receivers. For the multi-path Rayleigh 
fading channel considered here the interference statis- 
tics is strictly speaking high-degree %-squarc for which 
Gaussian statistics is a reasonable approximation for 
our purpose of deriving a rough estimate of the trade- 
off between the stability and information rate. 

According to the preceding analysis SIR ~ 
NB/(MC) in the case of BBFS channels and SIR ~ 
N/M in the case of NBFN channels, independent 
of the power constraint. Let us set the covariance 
matrix of the receiver noises to be vIm> the M-dim 
identity. Suppose the average transmission power is 
constrained to P is constrained and all the trans- 
mit and receive antennas are identical. We have 
SNR = P/(vMG). In the case of BBFS channels, 
the SINR is roughly 



SINR 



MC 
~~NB 



isMCy 1 
P ) 



and, in the case of NBFN channels, it is 



SINR 



M 



vMC 



P ) 



(7) 



(8) 



The optimal SINR is obtained by setting N ~ 
P(i/max (B, C)) -1 which is roughly P/(vB) since 
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FIG. 2: Single-layer pinholes 



C < 2B. The information rate R, given roughly by 
2~ 1 MCln (1 + SINR), achieves the maximum roughly 
equal to P/v (in nat) at MC > TVi? (for BBFS) or 
N > M (for NBFN). That is, the channel capacity is 
linearly proportional to the power and is achieved at 
the expense of statistical stability of signals. 

Consider the thermal noise power v = ksTB where 
fee is the Boltzmann's constant and T the temper- 
ature. Then the above result implies that the en- 
ergy cost per nat is P/R ~ fc^T (with a constant 
close to unity) which is an extension to the TR-MIMO 
Rayleigh fading channel of the classical result derived 
for a simple SISO channel (see also pl|). 



V. PINHOLE EFFECT 

Pinholes are degenerate channels that can occur in a 
wide family of channels, outdoor as well as indoor, see 
Figs. 2 and 3. While preserving the co-channel decor- 
relation, pinholes have has been shown to severely 
limit the degrees of freedom and reduce the channel 
capacity 0, p|. In this section, I present an 
analysis of the pinhole effect on TR in Rayleigh fad- 
ing channels to demonstrate similar effects on stability 
and information rate. 

First, let us consider the simplest case of single-layer 
pinholes as illustrated in Fig. 2. Let hW(w) be the 
N x K transfer matrix from the TRA to the pinholes 
and h^ 2 ) (lu) the K x M transfer matrix from the pin- 
hole to the M receivers at frequency lu. The combined 
channel can be described by H(w) = h^ 2 ) (Lu)h^ (a;) = 

[h^ (co)h^j (lo)] in which and h\j are assumed to 
be i.i.d. C/V(0, <7i) and CiV(0, 02), respectively. 

At frequency lu, the mean signals received are given 
by EHHtm = -/Vi^oio^m and the variance of the 




FIG. 3: Multi-layer pinholes 
signal at frequency lu received at receiver i is given by 

KM N 

fe=l n— 1 j—1 

K KM N 

+E E i^ ) i 2 EE E i^ ) iVni 2 E E iCi 2E i< ) i 2 

k=l 1=1 n=l j=l 

K N 

+E( e i^ ) i 2 ) 2 ki 2 (E e i4- ) i 2 ei^ ) i 2 ) 
k=i j=i 

re KN(MN + MK + l)cr 2 cr||/x| 2 . 

Taking into account the temporal aspect of the sig- 
nal as before we obtain the normalized variance of 
the signal to the leading order (N, K > 1) for BBFS 
channels 

VjttMCB- l (N- x + K- 1 ), Vj 

and for NBFN channels 

Vj&M(N- 1 + K- 1 ), Vj. 

The result indicates that there is an effective number 
of TRA-elements given by iV cff = NK(N + K)^ 1 , 
namely one half of the harmonic mean of N and K, 
so that Vj « MCB^N^ 1 for BBFS channels and 
Vj ps MN'fi for NBFN channels . The confirms the 
intuition that pinholes are choke-points that reduce 
the effective number of TRA-elements. 

The previous case without pinholes corresponds to 
the limiting case K — > 00. For a fixed K, however, 
the previous benefit of having large number of TRA 
elements (N 3> 1) disappears. The multiple antennas 
in TRA are essentially screened out by the pinholes 
and the effective number of TRA-elements becomes 
K. 

The same analysis can be applied to channels with 
(n — 1) layers of pinholes such as illustrated in Fig. 
3. Let Kk, k — 1, ..n — 1 be the number of fc-th layer 
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FIG. 4: Example of a simple graph 

pinholes. Let h( fe ) be the transfer matrix for the k- 
th stage channel whose entries are i.i.d. zero-mean, 
variance-Ofc Gaussian r.v.s. and let the transfer matri- 
ces of different stages be independent of one another. 

I now show that to the leading order 
(N, K\, .., K n _i ^> 2") the normalized variance 
of the signal is given by V « MCB~ 1 N~g where the 
effective number of TRA-element N e g is given by 

-l -l 

namely n^ 1 times the harmonic mean of 
N,Ki,--- ,K n -\. The effective number of pin- 
holes N p is always larger than N e s and this again 
confirms the limiting nature of the pinholes. I sketch 
the proof below. 

The calculation of the mean is straightforward: 
EHH'm = NK\---K n ^\. Let us analyze the sec- 
ond moment of entry a 

E^HE^m) (HHWV 

- E V h (n) h (n - 1] •••h (2) h (1) h {1) *h {2) * 

— m, ^ ' t M„' t i n ,i„_i " , i 3 ,h"'i2M"'h,h"'h,h 

y ... /,(«)* ™. \^ h (n-l)* 

*1 »*■**«, 

x..-hyr,hy?,hV\,h™.-.h<r? hV . mi 

? 3 ,t 2 i 2 ' ( l ^2' 2 1 ^3'-?2 Jn-Jn — 1 J n +l^ n 

I claim that according to the Gaussian rule the lead- 
ing order terms of the second moment correspond 
to the simple graphs in which the arcs, connecting 
(un)primed indices to (un)primed indices, are nested 
and are bound by the ladder edges, connecting un- 
primed indices to primed indices. This includes the 
graph corresponding to |EHH^m| which has no lad- 
der edges. A simple graph is illustrated in Fig. 4. 

The observation is proved by induction. When a 
new layer of pinholes, described by hS n+1 \ is added, 
the number of graphs is doubled: one half of them con- 
tain the ladder edges connecting h^ +1 ^ to f&,* and 



hj 1 * to h^" while the other half contain 

the arcs connecting to hk 1 4 and 

to ./ . A moment of reflection reveals that the 

new pair of arcs impose one more constraint on the 
ranges of the indices than the new pair of ladder edges 
and yield a linear factor of K n+ i or M to the n-th or- 
der graphs while the new pair of ladder edges yield 
a quadratic factor MK n +\ to the graph correspond- 
ing to |EHH^m| and K 2 +1 to the rest of the n-th 
order simple graphs. Note that the structure of sim- 
ple graphs is not changed by the new ladder edges 
and the total number of n-th order simple graphs is 
exactly n + 1 . 

Collecting the terms corresponding to the simple 
graphs we have 

n n— 1 
i=l 3=1 

71—1 71—1 

x(Y[K k + N^K 1 ---K i ---K n _ 1 ) 

fe=i i=i 

where Ki means that Ki is absent in the product. 
Dividing it by N 2 an( ^ accounting for the 

temporal aspect of transmission we obtain the claimed 
result. 

To calculate the SINR we can substitute N e g for 
N in J7J). Even with unlimited power supply one 
would do best by choosing the optimal power level 
P ~ BvN C K] higher power level is wasteful. The max- 
imal information rate is roughly BN c g which is still 
bounded by BN p . 

VI. CONCLUSIONS 

I have analyzed the time reversal of propagation in 
multi-path Rayleigh-fading MIMO-channels with or 
without pinholes. The focus of the analysis is the sta- 
bility condition, the multiplexing gain and their trade- 
off. The main results can be summarized as (i) that 
the stability holds when MC < N eS B for BBFS chan- 
nels and M <C N e g for NBFN channels where N e g is 
the effective number of TRA-elements and (ii) that 
the maximal information rate BN c g is attained when 
the power supply is set to the noise level times BN c g 
and when the stability condition is reversed. They 
are a significant extension of the previously discussed 
conditions for stability Q, 0, 0, 

I have computed the effective number of TRA- 
element iVeff under the condition that the numbers 
of TRA-elements and the pinholes of each layer are 
much greater than 2", with n — 1 being the number 
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of layers. In this extreme case the effective number 
of TRA-elements is asymptotically equal to n^ 1 times 
the harmonic mean of TRA-elements and the numbers 
of pinholes at all n — 1 layers. In the opposite case the 
graph analysis shows that the normalized variance of 
signal grows exponentially with the number of pin- 
hole layers and consequently the rate of information 
transmission diminishes exponentially fast. In other 
words, a long chain of independently fluctuating me- 
dia separated by a series of screens of sparse pinholes 
is detrimental to time reversal communication. 

The estimate NB for the TR information rate in the 
absence of pinh oles should be contrasted with the find- 
ing in , H3 , H2 , H3 that the capacity with CSI at 



the receiver with M receive antennas (but not at the 
N transmit antennas) scales like B min (M, N) In SN R 
at high SNR. Note that their result does not include 
the interference due to non-cooperating multiuser re- 
ceivers as is done here and is for narrow-band signals 
with C = 2B. 
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